Chandra Observations of WR 147 Reveal a Double X-ray Source 

Svetozar A. Zhekov^''^ and Sangwook Park^ 



ABSTRACT 

We report the first results from deep X-ray observations of the Wolf-Rayet 
binary system WR 147 with the Chandra HETG. Analysis of the zeroth order 
data reveals that WR 147 is a double X-ray source. The northern counterpart is 
hkely associated with the colhding wind region, while the southern component 
is certainly identified with the WN star in this massive binary. The latter is the 
source of high energy X-rays (including the Fe Ka complex at 6.67 keV) whose 
production mechanism is yet unclear. For the first time, X-rays are observed 
directly from a WR star in a binary system. 

Subject headings: stars: individual (WR 147 ) — stars: Wolf-Rayet — X-rays: 
stars — shock waves 

To Appear in The Astrophysical Journal Letters 



1. Introduction 



The Wolf-Rayet (WR) star WR 147 is a massive binary system that is bright in X-rays 
and is a composite radio source. Being the second closest WR star known allowed its emis- 
sion to be spatially resolved in the radio, near infrared (NIR) and optical. High-resolution 
radio observations showed that its southern component, WR 147S (the WN8 component in 
the bin ary), is a thermal source while its northern counterpart, WR 147N , is a non-thermal 



source (Abbott et al.lll986 



Williams et al.l 



1997 



Moran et al 



Skinner et al 



1989; Churchwell et al. 1992; Contreras et al. 1996 



19991 1. High- resolution NIR dat a revealed two source s 



separated by ^ 0'.'64 and WR 147N was classified as a B0.5V star (IWilliams et al.l 119971 ). 
An earlier spectral class, 08-09 V-111, was suggested f rom Hubble Space Te lescope observa- 
tions which spatially resolved the binary components (INiemela et al.lll998l ). At a distance 
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of 630 ± 70 pc to WR 147 (IChurchwell et al.lll992l ). the projected (or minimum) binary sep- 
aration is 403 ± 13 au. Since the binary components are massive stars, the interaction of 
their strong winds should fo rm a colhding-stellar wind (CSW) region (jPrilutskii fc Usov 



19761: 



Cherepashchuk I Il976l ). Interestingly, from the analysis of radio and IR images. 



Williams et al.l (jl997l ) concluded that the WR 147N radio peak is displaced by 0'.'07 south- 
ward from its IR counterpart. This value is within the uncertainties of the data but if it is 
real, the non-thermal radio source can be associated with the CSW region located in vicin- 
ity of the OB companion. The CSW scenario is supported by the analysis of the Chandra 
High Resolution Camera (HRC) image (although with a very limited photon statistics and 
having no spectral resolution) which concluded that the X-ray emission is spatially extended 
and it peaks north of the WN8 star altho ugh a deeper X-ray image is needed to accurately 
determine the degree of spatial extension (jPittard et al.ll2002l ). It is worth noting that con- 
tribution from the massive st ars in the binary cann ot be excluded since these objects are 



X-ray sources themselves (e.g. ICiidel fc Naze II2009I and the references therein). 



X-ray observations are very important for studying CSWs. They provide us with direct 
information about the physical conditions in the hot plasma behind strong shocks, which 
are likely the place where non-thermal radio emission forms. The earlier X-ray observations 
of WR 147 have revealed the presence of thermal emission from high temperature plasma 
(ICaillault et al.lll985l : ISkinner et al.lll999l ). but the CSW plasma characterist ics were more 



tight ly constrained only from the higher signal-to-noise XMM-Newton data (ISkinner et al. 
20071 ). Skinner et al. reported a detection of the Fe Ka complex at 6.67 keV and their 



analysis showed that the plasma temper ature is higher than the maximum temperature in 
the colliding-wind shock. Consequently, IZhekovl (120071 ) modeled the XMM-Newton spectra 
successfully in the framework of the CSW scenario but only at the requirement the stellar 
wind velocities were by a factor of 1.4 — 1.6 higher than their currently accepted values. 

The fact that WR 147 is a wide binary at a relatively small distance gives us an oppor- 
tunity to examine the CSW phenomenon in more detail that is to carry on spatially resolved 
X-ray studies. This motivated our deep observations with the Chandra High Energy Trans- 
mission Gratings (HETG). In this Letter we report the first part of our analysis (based on 
the zeroth order data) which reveals a double X-ray source in WR 147 . 



2. Observations and Data Reduction 

WR 147 was observed with Chandra in the period Mar 28 - Apr 10, 2009. The observa- 
tions in configuration HETG-ACIS-S were carried out in eight occasions {Chandra Obslds: 
9941, 9942, 10675, 10676, 10677, 10678, 10893 and 10897) that provided 7,300 zeroth order 



- 3 - 



counts in a total effective exposure of 286 ksec. The CIAO 4.1.2111 data analysis software 
was used through this study in combination with the Chandra calibration database CALDB 
v.4.1.3. Following the CIAO Science Threads, the original or combined event files were 
used for determining the source coordinates, source variability and spectral extractions. By 
default, the pixel randomization is switched off in grating data. 



3. Analysis and Results 

To have a full advantage of the total source counts in our analysis of the zeroth order 
data, it is necessary to combine all data sets. As a prerequisite to this process, it is important 
to estimate the accuracy of the source coordinates in the resultant total image. This was 
done in two steps. 

First, the CIAO command wavdetect was run on each of the eight observations and 
an average position of the source was calculated. The errors of the mean were propagated 
from the individual errors derived from wavdetect. Second, all the zeroth order data were 
merged eight times (CIAO command merge-all) with a different reference image. The source 
coordinates were determined for each merged image by running wavdetect. The maximum 
difference in the source coordinates from one merged image to another was 0!002 and (y.'02, 
correspondingly for the right ascension and declination. The source positions are given in 
Table [H We note that the center of the source in the merged data is practically the same 
as the one determined from averaging coordinates of the source centers derived for each 
individual observation. All this makes us confident that our analysis is not influenced by 
the choice of the reference image. The derived X-ray coordinates of WR 147 are in good 
agreement (within the la of the Chandra positional accuracy @) with the ones determined 
for this object in the optical [HST Guide Star Catalog) and in the radio (Tabled]). 



3.1. Images 

The considerably better photon statistics (7,300 cts) of the total zeroth order HETG 



image compared to that of the previous Chandra HRC image (~ 148 cts, iPittard et al.ll2002l ) 
gave us an opportunity to study in detail the spatial extent of the X-ray emission from this 
object. It also allowed us to carry on image analysis in different photon energy bands. 



^Chandra Interactive Analysis of Observations (CIAO), 
^http:/ /acs. harvard.edu/proposer/POG: §5.4.1 
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The relatively small spatial scale of the WR 147 emission in the radio and optical 
(C'ST — (y.'64 which corresponds to ~ 1.2 — 1.3 ACTS pixels) prompted the use of subpixel im- 
ages and image deconvolution techniques. We applied a de convolution procedure, based on a 



Lucylll974l ) , that was successf ully used in the 



maximum likelihood method (lRichardsonlll972 
imag e analysis of SNR 1987A, an o bject with a sniall sp atial extent of 



I'.'G ( Burrows et al. 



2000l : IPark et al.ll200l 2004, 2006; iRacusin et all 120091 1. The main result is that the mor- 
phology of the X-ray emission of WR 147 changes in different energy bands: its northern 
part is brightest in the (1.0-2.0 keV) range and gradually weakens at higher energies. This 
likely indicates that a soft X-ray source is located there while most of the hard energy pho- 
tons originate in the southern part of WR 147 . This pattern is clearly seen from the raw 
subpixel X-ray images as well as from the deconvolved ones (Fig. [1]) and is also illustrated by 
the sources relative brightness derived from direct photon counting (see lower right panel in 
Fig. [3]). It is thus conclusive that WR 147 is a double X-ray source, consisting of a northern 
(WR 147N ) and a southern (WR 147S ) part with different spectral characteristics. A sep- 
aration of ~ 0'.'60 between the positions of maximum emission of WR 147N and WR 147S is 
directly measured from the deconvolved (1.0 - 2.0 keV) image. 

The existence of a spatial separation between WR 147N and WR 147S gets additional 
support from measuring the WR 147 coordinates in different energy bands. We ran the CIAO 
command wavdetect on the eight merged images. We see that there is an offset of ~ O'.'Sl 
between the source coordinates in the soft (1.0 - 2.0 keV) and hard (6.0 - 8.0 keV) energy 
bands (Tabled]). Interestingly, the source coordinates in the latter almost perfectly coincide 
with the ones in the optical and radio. It is thus conclusive that the X-ray hard-photon 
emitter (WR 147S ) is in fact the WR star itself (a WN8 object). Furthermore, a direct 
comparison of the WR 147 images in the radical and X-rays reveals a nice correspondence 
between the locations of the radio and X-ray WR 147N sources (Fig. [1]) which is suggestive 
that this is the CSW region in this binary system. 



3.2. Undispersed Spectra 



Guided by the results from the image analysis, we extracted the X-ray spectra of 
WR 147 , WR 147N and WR 147S from the total (merged) zeroth order data by using 
the CIAO script specextract (for extraction regions see Fig. [21 ). We fitted th e spectra with 
standard as well as custom models in version 11.3.2 of XSPEC (lArnaud 1119961 ). We note that 



•^The radio data are from the NRAO Archive. The National Radio Astronomy Observatory is a facihty 
of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. 
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due to the small spatial separation between WR 147N and WR 147S their spectra are not 
completely disentangled from each other. In our upcoming analysis of the first-order data, 
we will consider a more elaborated modeling of the undispersed spectra in conjunction with 
the dispersed ones. We adopted a simpler approach here which is nevertheless instructive 
for the emerging physical picture. 

As seen fro m Fig. [3l the total s pectrum of WR 147 is very similar to that obtained with 
XMM- Newton (jSkinner et al.l 120071 ): it is heavily absorbed, having most of its X-ray counts 
at energies E > 1 ke V. For consistency with the previou s studies, we adopted the same set of 



WN abundances (see lSkinner et al.ll2007l : IZhekovll2007l ). A simple IT shock model [vpshock 



in XSPEC) gave a similar quality of the fit as for the XMM-Newton spectra of WR 147 . The 
X-ray absorption and the plasma temperature (with 90% confidence intervals in brackets) 
were = 2.5[2.3-2.9] x 10^^ cm'^, kT = 2.9[2.3-3.2] keV, and the de rived abundarices fe ll 
within the 90% confidence intervals of the XMM-Newton spectral fits (jSkinner et al.l 120071 ). 
The total observed flux was Fx(0.5 — 10 keV) = 1.28 x 10^^^ erg s~^, about 15% smaller 
than that from the XMM-Newton spectra. A better constraint on this difference would be 
possible after completing upcoming analysis of the first-order HETG data. 

Exploring the physical picture that identifies WR 147N with the CSW region in the 
binary system, the X-ray spectra of WR 147N and WR 147S were fitted simultaneously and 
they shared the same WN abundances set. The latter is justified by that the shocked WN 
wind dominates the X-ray emission of the CSW region (jZhekovl 120071 : for the use of discrete 
temperature models see § 5.2 therein). The IT shock model with individual postshock 
temperatures for WR 147N and WR 147S and a common X-ray absorption gave a poor 
quahty of the fit (reduced ~ 1-5)- An acceptable fit was obtained if the spectra had 
individual X-ray absorption (see Fig. [3] and Table [2]). It should be emphasized that in both 
sources the plasma temperature well exceeds 1 keV. 

As in the image analysis, we see that the northern source is softer than its southern 
counterpart. Also, the latter (the WN8 star) is the place where the Fe Ka complex at 
6.67 keV comes from (see the inset in Fig. [3]). W R 147S is s u bject to an excess X-ray 
absorption compared to WR 147N . Applying the iGorensteinI (119751 ) conversion (NH = 
2.22 X 10^^ Ay cm~^), the absorption towards WR 147N almos t perfectly correspon ds to the 
optical extinction of WR 147 (A^ = 11.6 mag. Ay = A^/1.11; Ivan der HuchtlboOlh . 



We thus note that the results from image and spectral analysis seem consistent with a 
physical picture where WR 147N resides in the CSW region of the binary, while the X-ray 
source WR 147S is likely located deeper in the WR wind. 
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3.3. MARX Simulations 

We used version 4.4 of the MARX0 software to simulate the observational situation with 
WR 147 . We ran a 286-ksec HETG exposure for a X-ray source composed of two point 
sources with spectral characteristics corresponding to those of WR 147N and WR 147S (Ta- 
ble The sources were located at the positions of maximum emission for the northern and 
southern sources as in the (1.0-2.0 keV) deconvolved image of WR 147 . We note a nice 
correspondence between the simulated and observed data (Fig. [1]). 

All this makes us more confident about the results from our analysis, that illustrate the 
superior capabilities of the Chandra observatory even when working on their very edge. 



4. Discussion 



Analysis of the zeroth order HETG data showed that Chandra resolved WR 147 into a 
double X-ray source. Its two counterparts, WR 147N and WR 147S , are most likely identified 
correspondingly with the CSW region and the WN8 star in this wide binary system. We 
recall that the maximum temperature in the CSW region is determined by the terminal wind 
velocities in the binary and for WR 147 their currently accepted values cou ld not provide 
the plasma temperature deduced from the X-ray spectra (jSkinner et al.l 120071 ). This is why, 
successful CSW models of WR 147 were possible only for wind velocities being a factor of 
1.4 — 1.6 higher (jZhekovl 120071 ). We note that all this is based on analysis of the unresolved 
XMM-Newton data which assumed that the CSW region is responsible for the total X-ray 
emission from WR 147 . The high resolution Chandra data resolve this discrepancy by 
revealing two sources with different spectral characteristics. 

We used the CSW model with the nominal stellar wind parameters {Vwr = 950 km s~^, 
MwR = 4x10-5 Mq yr-i; Vq = 1 600 km s-\ M p = 6.6x10"^ M© yr'^; [MoVq/ MwrVwb] = 



0.028; for CSW model details see IZhekovll2007f ) to fit the X-ray spectrum of WR 147N . As 
seen from Fig. [3], the shape of the observed spectrum is perfectly matched, that is there 
is NO temperature discrepancy between the model and observations any more. But, the 
model overestimates the total luminosity of WR 147N by a factor of ~ 16. We note that 
this discrepancy is in general trackable. The CSW X-ray luminosity scale s with the mass 
loss rate M, wind velocity v and binary separation as Lx oc M'^v~^D~^ (ILuo et al.l Il990l : 
Myasnikov fc Zhekovl 119931 ). If the mass- loss rates are factor of ~ 2 or even more lower 



than assumed because of their intrinsic uncertainties or because the winds are clumped (e.g. 



http://space.mit.edu/CXC/MARX/ 
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Crowtherl 120071 ) . this can account for most of the luminosity mismatch, and the unknown 
orbital inclination (larger binary separation) would account for the difference. All this will 
be explored in detail in our upcoming combined modeling of the first-order spectra and the 
undisperesed ones. 

While this gives more confidence that WR 147N resides in the CSW region of the binary 
system, an identification of WR 147N with the OB star in the system cannot be completely 
ruled out. Our analysis of the X-ray spectral lines (from the first-order HETG spectra) will 
be crucial in this respect. For example, if suppressed forbidden line in the He-like triplets is 
found (an indication of high density or strong UV emission), it will favor the OB star origin 
for WR 147N . Unfortunately, the luminosity of the OB companion is poorly constrained 



(ILepine et al.l l200ll ) which prevents estimating its contribution to the X-ray emission of 
WR 147N based on the Lx - l-iboi relation. But, the OB stars are in general soft X-ray 
sources with plasma temperatures less than 1 keV (see §4.1.2 and §4.3 in the review paper of 
Giidel fc Naze II2009I ). Given the high X-ray absorption to WR 147 (§ 13. 2p . we thus anticipate 
that this contribution will not be appreciable un less the OB star in WR 147 is a rare hot 
magnetic object (see §4.6 in loiidel fc Naze Ibood ). 

On the other hand, the fact that the WN8 star is a hard-energy X-ray source (WR 147S ) 
is indeed puzzling. Skinner et al. (2002a,b) have discussed in some detail possible mech- 
anisms for X-ray production in WN stars when the high-energy tail was established with 
certainty in the X-ray emission from presumably single WR stars: WR 6 and WR 110. 
These mechanisms include: instability-driven wind shocks (IDWS); magnetically confined 
wind shocks (MOWS); wind accretion shocks; colliding wind shocks (including the case of 
stellar wind shocking onto a close companion); non-thermal X-ray emission. We note that 
none of these mechanisms finds solid observational support for the moment and each of 
them has its own limitations and caveats. For example, the IDWS should be considerably 
soft X-ray emitters while the MOWS might be able to provide the high plasma temperature 
observed in WR 147S . But for the massive stellar winds in WRs, this requires a relatively 
strong global magnetic field in the star itself. At present, such fields have not been reported 
for Wolf-Rayet stars. Presence of a strong magnetic field may also suggest fiare activity but 
no short-term variability is detected in any of the eight HETG observations (the CIAO tool 
glvary retur ns a variability index of zero). The same was found in earlier X-ray observations 
of WR 147 (jSkinner et al.lll999l : ISkinner et al.l 120071 ). Presence of a compact (or a normal 
star) companion is intriguing but this may suggest that we have quite a rare opportunity 
to observe a triple system that initially consisted of at least two or even three massive stars 
and in the latter case the most massive one has already exploded. 



Finally, although the mechanism responsible for the high-energy X-ray production is 
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uncle ar it may not be uni que to WR 147 . Apart from earlier detections (WR 6 and WR 
110), ISkinner et al.l (120091 ) report a hard-energy tail in the X-ray spectra of several more, 
presumably single, WN stars. It is worth noting that all these WNs have a subtype different 
from WR 147 while a star of its subtype, WR 40 (a WN8 star), was not detected with 
XMM-Newton (iGosset et al.ll2005l ). Hopefully, future observations will facihtate solving the 
mystery of this puzzling phenomenon: how common the high-energy tail is for single WR 
stars and how it correlates with their subtype. But, WR 147 is unique in the sense that 
for the first time we detect X-rays directly from a WR star in a binary system. Moreover, 
the upcoming analysis of the first-order HETG spectra will give us a chance to measure 
line parameters (widths, centroids, fluxes) in the X-ray spectrum intrinsic to a WR star 
(WR 147S ), thus, to study its X-ray plasma characteristics in detail. 



This work was supported by NASA through Chandra grant GO9-0013A to the University 
of Colorado at Boulder, and through grant G09-0013B to the Pennsylvania State University. 

Facilities: CXO (HETG, ACTS). 
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Fig. 1. — Examples of WR 147 X-ray images (pixel size of 0'.'123) in a linear scale by 
rows: (1) raw images; (2) deconvolved ones; (3) images from MARX simulations and a radio 
image. RA(J2000) and Dec(J2000) are on horizontal and vertical axes, respectively. The 
optical positi on of WR 147 (HST GSC) is m arked by a circle and the position of the southern 
radio source (IContreras fc Rodrigueg|l999l ) is marked by a square. The scale of each image 
is illustrated by the (1" x 1") square in its lower right corner. The 3.6-cm radio image of 
WR 147 is from a VLA observation on June 28, 1999. The bright southern source is the WN8 
star while the fainter emission to its north is the non-thermal source (presumably arising in 
CSWs). Overlaid are the contours (linearly spaced by 0.1 of the maximum emission) of the 
deconvolved X-ray image in the (1.0-2.0 keV) energy range. The radio and X-ray images are 
in relative units normalized to their corresponding maximum brightness. They were aligned 
that the brightness peak of the southern X-ray source coincided with that of the southern 
radio source. The region boxes (in green) were used for estimating the relative brightness 
of WR 147N and WR 147S (see Fig. [3]). Their border line defines how the two sources were 
separated for the spectral extractions (Fig. [2]). 
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Fig. 2.- 



The extraction regions for tlie total WR 147 , WR 147N and WR 147S spectra. 
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Fig. 3. — WR 147 background-subtracted spectra (rebinned to have a minimum of 20 
counts per bin). Upper left panel: The total spectrum overlaid with the best-fit IT shock 
model ix^/dof = 269/230). Prominent emission lines are marked. Upper right panel: the 
WR 147N and WR 147S (in red) spectra and the IT shock models with individual X-ray 
absorption. The inset shows the Fe complex at 6.67 keV. Lower left panel: the X-ray 
spectrum of WR 147N overlaid with the CSW model {x^/dof = 98/80). Lower right panel: 
The WR 147N -to-WR 147S relative brightness defined in different energy ranges (marked 
with horizontal bars) by counting the number of photons in the two region boxes as shown 
in Fig. [H 
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Table 1. WR 147 Positions: X-ray, optical, radio 
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Note. — The source position by rows: (1) the mean 
from the eight Chandra observations; (2) from the 
merged images; (3, 4) from 'filtered' merged images; 
(5) optical; (6, 7) radio. 



T he radio coord i nates (jContreras fc Rodrigue2 



19991 : IWilhams et al.l Il997f ) are for WR 147S (the 
WN8 star) converted from B1950 into J2000 
with the HEASARC coordinate converter 



( http: / /heasarc. gsfc.nasa.gov/cgi-bin/Tools/convcoord/convcoord. pi ) . 
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Table 2. IT Shock Model Results 



WR 147N WR 147S 



xVdof 309/256 

Nh (1022 cni-2) 2.28 [2.08 - 2.57] 3.83 [3.51 - 4.20] 

kT (keV) 1.78 [1.52 - 1.98] 2.36 [2.12 - 2.56] 

Ne 19.7 [0.0 - 76.2] 

Mg 3.3 [0.5 - 6.7] 

Si 5.1 [3.9 - 7.2] 

S 6.7 [5.5 - 7.8] 

Ar 8.4 [5.8 - 9.8] 

Ca 8.3 [4.3 - 12.4] 

Fe 10.0 [7.7- 11.8] 

Fx'' (10-12 ergs cm-2 s^i) 0.324 (4.4) 0.921 (12.3) 



Note. 



Brackets enclose 90% confidence intervals. All abun- 



dances are with respect to their solar values (lAnders fc Grevesse 
19891 ). The fixed in the fit abundances are: H= 1, He= 25.6, 
C= . 9, N= 140, = 0.9, and Ni= 1 (for details see lskinner eTaL 



20071 : IZhekovll2007f ) 



"^The observed X-ray fiux (0.5 - 10 keV) followed in parentheses 
by the unabsorbed value. 



